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Graphene, a one-atom thick zero gap semicon- 
ductor [l|, |2|], has been attracting an increasing 
interest due to its remarkable physical properties 
ranging from an electron spectrum resemb ling 
relativistic dynamics [3j, 3, [g, [g|, 0, [g|, [9, KLOL IllL LL2I ] 
to ballistic transport under ambient conditions 
d 0? S, 3] . The latter makes graphene a promis- 
ing material for future electronics and the re- 
cently demonstrated possibility of chemical dop- 
ing without significant change in mobility has im- 
proved graphene's prospects further [13]. How- 
ever, to find optimal dopants and, more generally, 
to progress towards graphene-based electronics 
requires understanding the physical mechanism 
behind the chemical doping, which has been lack- 
ing so far. Here, we present the first joint experi- 
mental and theoretical investigation of adsorbates 
on graphene. We elucidate a general relation be- 
tween the doping strength and whether or not 
adsorbates have a magnetic moment: The para- 
magnetic single NO2 molecule is found to be a 
strong acceptor, whereas its diamagnetic dimer 
N2O4 causes only weak doping. This effect is re- 
lated to the peculiar density of states of graphene, 
which provides an ideal situation for model stud- 
ies of doping effects in semiconductors. Further- 
more, we explain recent results on its "chemical 
sensor" properties, in particular, the possibility 
to detect a single NO2 molecule [13]. 

Controlling the type and the concentration of charge 
carriers is at the heart of modern electronics: It is the 
ability of combining gate voltages and impurities for lo- 
cally changing the density of electrons or holes that al- 
lows for the variety of nowadays available semiconduc- 
tor based devices. However, the conventional Si-based 
electronics is expected to encounter fundamental limita- 
tions at the spatial scale below 10 nm, according to the 
semiconductor industry roadmap, and this calls for novel 
materials that might substitute or complement Si. Being 
only one atomic layer thick, graphene exhibits ballistic 
transport on a submicron scale and can be doped heav- 
ily — either by gate voltages or molecular adsorbates — 
without significant loss of mobility [1]. In addition, later 
experiments [Hj demonstrated its potential for solid state 



gas sensors and even the possibility of single molecule de- 
tection. 

We show, that in graphene, aside from the donor- 
acceptor distinction, there are in general two different 
classes of dopants — paramagnetic and nonmagnetic. In 
contrast to ordinary semiconductors, the latter type of 
impurities act generally as rather weak dopants, whereas 
the paramagnetic impurities cause strong doping: Due to 
the linearly vanishing, electron-hole symmetric density of 
states (DOS) near the Dirac point of graphene, localised 
impurity states without spin polarisation are pinned to 
the centre of the pseudogap. Thus impurity states in 
graphene distinguish strongly from their counterparts in 
usual semiconductors, where the DOS in the valence and 
conduction bands are very different and impurity levels 
lie generally far away from the middle of the gap. 

Impurity effects on the electronic structure of two- 
dimensional systems with Dirac spectrum were investi- 
gated in detail in connection with the problem of high- 
temperature superconductivity in copper-oxide com- 
pounds, and STM visualization of the order parameter 
around impurities is one of the most straightforward evi- 
dences of the (i-wave pairing in these systems [14| . For the 
case of a strong enough impurity potential, the formation 
of quasilocalized electron states in the pseudogap around 
the Dirac point is expected [14| , whereas a weak potential 
will not lead to the formation of quasilocalized states at 
all. This means that, in general, in this situation one can 
hardly expect a strong doping effect which requires ex- 
istence of well-defined donor (or acceptor) levels several 
tenths of electron volt away from the Fermi level. How- 
ever, if the impurity has a local magnetic moment its en- 
ergy levels split more or less symmetrically by the Hund 
exchange, of the order of 1 eV, which provides exactly 
the situation favorable for a strong doping. To check 
this assumption we have chosen the NO2 system form- 
ing both paramagnetic single molecules and diamagnetic 
dimers N2O4. Note that magnetism of impurity atoms in 
graphene is interesting not only due to its effect on the 
electronic properties but also itself. As was shown re- 
cently [lH ferromagnetism based on spin polarization of 
sp-electrons in narrow impurity bands should be charac- 
terized by much higher Curie temperatures than typical 
for conventional magnetic semiconductors. 



2 



Doping and gas sensing effects in graphene related sys- 
tems attracted a lot of research activity in the last years. 
Already the first experiments with graphene showed the 
possibility of inducing charge carriers to this material 
by the adsorption of various gases including NH3, H2O 
and NO2 [1]. Hall effect measurements proved, that NH3 
induces electrons, whereas the latter two types of ad- 
sorbates result in holes as charge carriers. Furthermore, 
they showed that the detection of those gases at remark- 
ably low concentrations or for NO2 even in the extreme 
limit of single molecules is possible [13j . 

Carbon nanotubes (CNT) being rolled up graphene 
sheets exhibit similar doping effects on gas exposure [l6[ 
and stimulated first principles studies of these systems: 
DFT calculations for NO2, H2O and NH3 on nanotubes 
revealed possible physisorbed geometries on nondefective 
CNTs and developed a "standard model" to interpret 
this doping [13, H, EH: By considering of Mulliken or 
Lowdin charges of the adsorbed molecules, NO2 is found 
to accept 0.1 e~ per molecule from the tube, whereas one 
NH3 molecule is predicted to donate between 0.03 to 0.04 
e~ [HI, El. However, this "standard model" for CNTs 
fails for graphene, especially in explaining the qualitative 
difference between NO 2 and the other adsorbates. 

In this work we study the NO2 and N2O4 adsorbate ef- 
fects discussed above by combining ab-initio theory with 
transport measurements. Theoretically, the electronic 
properties of carbon related materials and the predic- 
tion of adsorbate geometries from first principles can be 
well addressed by the density functional theory (DFT). 
Although van der Waals forces are ill represented in the 
local density approximation (LDA) as well as in gradient 
corrected exchange correlation functionals (GGA) result- 
ing in over- and underbonding, respectively [20], for polar 
molecules like NO 2 these errors are of minor importance: 
In similar studies of NO2 on nanotubes both function- 
als yield qualitatively the same predictions of adsorp- 
tion energies and geometries [21|. Here, we apply both 
functionals, so that we obtain upper and lower bounds 
for adsorption energies and related structural properties, 
such as the equilibrium distances between molecules and 
graphene plane. 

All DFT calculations were carried out with the 
Vieanna Ab Initio Simulation Package (VASP) [22[ using 
projector augmented waves (PAW) [23|,|24| for describing 
the ion cores. The corresponding plane wave expansions 
of the Kohn-Sham orbitals were cut off at 875 eV in the 
GGA [HHH and at 957 eV in the LDA-calculations. In 
this periodic scheme, single NO2 and N2O4 adsorbates 
are modelled in 3 x 3 and 4x4 graphene supercells, 
respectively. The ionic configurations presented in this 
letter are fully relaxed, i.e. all forces being less than 
0.02 eV A -1 , and the convergence of subsequent total en- 
ergy calculations is guaranteed by applying the tetrahe- 
dron method with Blochl corrections on T-centered k- 
meshes denser than 30 x 30 x 1, when folded back to the 
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Figure 1: Left: Spin-polarised DOS of the graphene super- 
cells with adsorbed NO2, a) - b),and DOS of graphene with 
N2O4, c) - e), in various adsorption geometries. The energy 
of the Dirac points is defined as Ed = 0. In the case of NO2 
the Fermi level Ef of the supercell is below the Dirac point, 
directly at the energy of the spin down POMO, whereas for 
N2O4 Ef is directly at the Dirac points. Right: Adsorption 
geometries obtained with GGA. The carbon atoms are printed 
in blue, nitrogen green and oxygen red. 



single graphene Brillouin zone. 

In the spirit of Ref. [27], the DOS obtained in our 
DFT calculations are the central quantities in the follow- 
ing discussion of the adsorbate effects on the electronic 
properties of the graphene sheets. In particular, the oc- 
currence of a strong acceptor level due to individual NO2 
molecules is predicted in this way, as well as experimen- 
tally confirmed. 

Gaseous NO2 stands in equilibrium with its dimer 
N2O4 giving rise to various different adsorption mech- 
anisms on graphene. For both, we obtained possible ad- 
sorption geometries as depicted in Fig. [1] right. The 
corresponding adsorption energies in GGA are 85meV 
(a), 67meV (b), 67meV (c), 50meV (d) and 44meV (e) 
per molecule with sheet-adsorbate distances of 3.4 — 3. 5A 
for the monomer and 3.8 — 3.9A for the dimer. As usual, 
LDA yields higher adsorption energies - approximately 
169 — 181 meV for the monomer and 112 — 280 meV for 
the dimer - and favours the adsorbates by 0.5 — 1 A nearer 
to the sheet. 

The spin-polarized DOS of the supercells containing 
NO2, shown in Fig. [1] a) and b), reveals a strong ac- 
ceptor level at 0.4 eV below the Dirac point due to these 
adsorbates in both adsorption geometries. 

The molecular orbitals of NO2 correspond to flat bands 
and manifest themselves as peaks in the DOS. The ener- 
gies of these peaks are virtually independent of the ad- 
sorbate orientation. Most important for doping effects is 
the partially occupied molecular orbital (POMO) of NO2, 
which is split by a Hund like exchange interaction: The 
spin- up component of this orbital is approximately 1.5 eV 
below the Dirac point and fully occupied, as it is also for 
the case of free NO 2 molecule. The spin down component 
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of the N0 2 POMO is unoccupied for free N0 2 , but 0.4 eV 
below the Dirac point in the adsorbed configuration . 
Hence, it can accept one electron from graphene. 

In contrast to the paramagnetic monomer, the dimer, 
N2O4, has no unpaired electrons and is diamagnetic: 
on formation from two monomers the two POMOs hy- 
bridize. The resulting bonding orbital is the highest occu- 
pied molecular orbital (HOMO), whereas the antibond- 
ing combination forms the lowest unoccupied molecular 
orbital (LUMO) for free N 2 4 . The possibility of dop- 
ing effects due to adsorbed dimers has been investigated 
using the DOS depicted in Fig. [T]c) - e). 

Again, the molecular orbitals of the adsorbates are rec- 
ognizable as sharp peaks in the supercell DOS. One finds 
that the HOMO is in all cases more than 3 eV below the 
Fermi level and therefore does not give rise to any dop- 
ing. However, the LUMO is always quite near to the 
Dirac point, i.e. between 1 meV and 66meV above it 

These initially empty orbitals can be populated by the 
graphene electrons due to thermal excitations and act 
consequently as acceptor levels. Thus both, N2O4 and 
NO2, give rise to p-type doping of graphene — with 
one decisive distinction: The affinity of the paramag- 
netic monomer to accept electrons from graphene is much 
stronger than for the dimer, manifesting itself in the en- 
ergies of the acceptor levels being well below and near 
the Dirac point, respectively. This confirms exactly our 
estimations that magnetism is crucial for a strong dop- 
ing effect. Without spin polarization, the quasilocalized 
states are typically formed in the pseudogap of the host 
DOS: The resonant energy of an impurity state is approx- 
imately located, where the real part of the denominator 
of the scattering T-matrix vanishes [3]. Here, however, 
it is Hund exchange splitting which creates acceptor (or 
donor) peaks at distances of order of half of an eV which 
is optimal for strong doping, as we observed for the case 
of N0 2 . 

Experimentally, the existence of these two distinct 
types of impurity levels was confirmed by combining elec- 
tric field effect and Hall measurements at different ad- 
sorbate concentrations. With the experimental setup de- 
scribed in [13[ we measured the Hall resistance R xy of an 
N02-exposed graphene sample as a function of the gate 
voltage Vq, where different adsorbate concentrations are 
achieved by repeated annealing of the sample at 410 K. 

At magnetic fields B = 1 T and room temperature 
T = 290 K we obtained l/R xy versus Vq as depicted in 
Fig. [2j These measurements exhibit two characteristic 
features. 

For gate voltages below — 10 V the curves are shifted 
in parallel towards lower 1/R xy values when increasing 
the absorbate concentration. Furthermore, the peak cor- 
responding to i? xy — > is shifted to higher gate volt- 
ages, broadens and the minimal |l/i? xy | values increase 
with more adsorbates on the sheet. We argue that this 
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Figure 2: The inverse Hall resistance 1/R xy is shown as a 
function of the gate voltage Vb for a graphene sample with 
different concentrations of NO 2 on top. The solid lines are the 
experimental results with the margenta curve corresponding 
to the highest and the green curve to the lowest concentration 
of adsorbates. The dotted lines are the simulations. 



behaviour, as a whole, is a direct evidence of the two 
distinct acceptor levels (See Fig. [3]). 

Depending on the adsorbate concentration and the 
gate voltage there are in general three different types 
of charge carriers contributing to the Hall effect: elec- 
trons and holes in graphene as well as electrons in the 
impurity states. Assuming an equal electron and hole 
mobility in graphene, the inverse Hall resistance is given 
by 1/Rxy = e ^t c n ^ n P T pr where n (p) is the density of 
electrons (holes) in graphene, c is the density of elec- 
trons in impurity states and \i their mobility in units of 
the graphene electron mobility. 

We determined \i « 0.08 from our electric field ef- 
fect measurements and the only adjustable is parame- 
ter density of impurity states N[ mp (E) yielding c via the 
Fermi distribution. As the gate voltage Vq = acr is di- 
rectly related to the total charge density of the sample 
o~ = e(c + n — p), where the prefactor a is determined by 
substrate properties as described in [l[, we can simulate 
the Hall resistance as a function of Vq by assuming an 
explicit form of the impurity DOS. 

It turns out, that for reasonable agreement of all ex- 
perimental curves with the simulations Ni mp (E) has to 
be peaked around two distinct energies, E\ < — 300 me V 
and E2 ~ — 60meV. Therefore, the experiment yields two 
types of acceptor levels, one rather close to the Dirac 
point and the other well below, in agreement with the 
DFT predictions. The simulations shown in Fig. [2] are 
according to these impurity level energies, in particu- 
lar N imp (E) = Cl S(E - E x ) + c 2 S(E - E 2 ) [31]. Fit- 
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Figure 3: The inset shows the DOS of graphene (magenta) 
with the strong acceptor level (blue) and the weak acceptor 
level (green) near the Dirac point. These impurity levels man- 
ifest themselves in the experimental l/i? xy versus Vb curves 
shown in the main panel. After long annealing (84h) one ob- 
tains the brown curve, whereas red curve corresponds to a 
higher adsorbate concentration on the sheet. The deep ac- 
ceptor level causes a solid shift at all Vb, while the acceptor 
level close to the Dirac point gives rise to an additional shift 
only at larger positive gate voltages. 
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Figure 4: The measured adsorbate concentrations are shown 
as a function of annealing time. Those were obtained by iden- 
tifying the NO2 (N2O4) concentrations with the density of 
impurity states ci (02). The latter have been determined by 
fitting the simulated 1/R xy versus Vb curves to our experi- 
ments. 



ting the number of impurity states c\ {pi) at — 300 me V 
(— 60meV) to each experimental curve yields the values 
for ci and C2 depicted in Fig. [H 

One sees that the adsorption of monomers and dimers 
is similarly likely, so that understanding the interplay of 
adsorbate and gate voltage effects on the transport prop- 
erties of NO 2 exposed graphene requires accounting for 
monomers as well as dimers. Note also, that the simi- 
lar adsorption and annealing rates for NO2 and N2O4 in 
experiment correspond nicely to the fact, that in both, 
GGA and LDA, the adsorption energies of the monomer 
are quite close to those of the dimer. 

In conclusion, we have elucidated the microscopic ori- 
gin of doping effects due to molecular adsorbates on 
graphene and work out the crucial relation between mag- 
netic moments and doping strength in this context. In 
particular, we prove the occurrence of a strong accep- 
tor level due to single NO2 molecules which is associated 
with the formation of local magnetic moment. Especially 
this strong acceptor level is the origin of many interest- 
ing phenomena and a promising candidate for tailoring 
the electronic and magnetic properties of future graphene 
devices. Besides explaining the effect of single molecule 
detection reported in [13[ it is capable of controlling the 
occupancy of flat impurity bands in graphene near the 
Dirac point and can give rise to exchange scattering. 
This earns future attention, as it can result in strongly 
spin-polarized impurity states [28] and guides a possible 



pathway [15] to high temperature magnetic order in this 
material. 
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